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~~i Abstract 

c3 . 

OnThe measurement of large scale patterns or anisotropies in the arrival direction of high energy cosmic rays is an important step 
^ fowards the understanding of their origin. Such measurements rely on an accurate estimation of the detector relative exposure in 
5^ each direction on the sky: the coverage map. To reach an accuracy on the determination of this map below the 1% level one must 
^ properly identify and correct for all the environmental effects that may induce variations in the detector exposure as a function 
d pf time. In an approach, similar to the one used in anti-sidereal time analysis, we propose a method to empirically estimate and 
""^correct for those effects under the hypothesis that seasonal and diurnal variations can be factorized. We tested this method using 
k model ground detector of cosmic ray air showers, whose aperture varies due to the dependence of the air shower development on 
the atmospheric conditions. 
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1. Introduction 



' Large scale anisotropy studies play an important role 
in the search for the origin of the highest energy cosmic 
rays. While some results [1] have favored a scenario with a 
~4% dipole pointing towards the Galactic Center and the 
Cygnus region at EeV (1 EcV = 10^* eV) energies, analyses 
with Auger data [2] have shown consistency with isotropy 
in this region of the sky. A series of large scale complemen- 
tary analyses with the first two years of Auger data was per- 
formed in [5] , focusing on right ascension modulations and 
carefully accounting for systematic effects, showing as well 
no indication of significant anisotropy at EeV energies. At 
the PeV (1 PcV = 10^^ eV) range, the Rayleigh formalism 
applied to right ascention modulations of the KASCADE 
experiment also shows no hints of anisotropies, both for 
light and heavy primaries [4] . 

In the 10^^ — 10^^ eV range, underground muon and Ex- 
tensive Air Shower (EAS) detectors were able to shed some 
light into the nature of large scale anisotropies at these en- 
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ergies |5l6l7l8j . For energies ^ eV, explored mainly by 
shallow underground muon detectors, cosmic rays are be- 
lieved to be modulated by solar magnetic fields, whereas 
for energies 10^"* eV, dominated by EAS arrays, the pri- 
maries will feel essentially the local structure of the galac- 
tic magnetic field. Amplitudes of right ascention modu- 
lations measured in these two energy ranges are ^ 10~* 
and ~ 10~^, respectively. An Earth based detector mov- 
ing through an isotropic cosmic ray plasma at rest will reg- 
ister an intensity modulation due to the Earth's rotation 
around its axis, the Compton-Getting (CG) effect [5]. For 
the movement of the Earth around the Sun, the CG effect 
should give rise to a solar diurnal modulation. In fact, solar 
diurnal variations measured at ^ 10 TeV [7] (1 TeV = 10^^ 
eV) are consistent with such an expectation. In a similar 
way, in the cosmic ray plasma rest frame, the movement of 
the whole solar system around the galactic center could give 
rise to a similar CG effect, however, now the modulations 
should be seen in sidereal time. Analysis in the ~ 300 TeV 
range [10] showed no such modulations, indicating that at 
theses energies, cosmic rays corotate with the solar system 
trapped by the local galactic field inside the arm. 

In the energy interval 10^* - 10^^ eV, linking the two 
ranges discussed in the first and second paragraphs of this 



section, the EAS-TOP collaboration has reported a solar 
time modulation of amplitude ~ 3 x lO""*, in very good 
agreement with the expected value from the CG effect pT| . 
Up to 10^"' eV, a clear signal of sidereal modulation at the 
~ 3 X 10~* level were also seen by this experiment, and 
above 300 TeV the amplitudes are not significant, showing 
consistency with the Tibet results. 

Although one-dimensional methods, such as the first har- 
monic analysis (in the case of a uniform right ascention ex- 
posure detector) of the Rayleigh formalism, the differential 
East- West method [8] and the Fourier transform method of 
modified times fl2| , are quite powerful in identifying distor- 
tions in the right ascention or sidereal time distributions, a 
complete caracterization of large scale anisotropics rely on 
the accurate knowledge of the detector relative exposure in 
all visible directions of the sky, the two-dimensional cover- 
age map. The proper estimation of this map is a very chal- 
lenging experimental task. All changes in the local (Earth 
based) running conditions of the detector that may impact 
on the exposure must be identified and taken into account. 
A perfect coverage map should contain all variations in- 
duced by local effects and none of those coming from the 
sky. In such an ideal case, the excess maps, i.e. the ratio of 
the observed signal over the estimated coverage in each di- 
rection of the sky, will only show the true sky anisotropics. 

The exposure in a given direction of the sky will be a 
function with a geometrical part depending on the local de- 
tection angles: the zenith angle 9 and the azimuth 0, and 
a part depending on time. Even though our knowledge on 
extensive air shower development through the atmosphere 
has clearly increased in the last decades, it is still incom- 
plete and Monte Carlo approaches, in addition to being 
computationally demanding, have a range of applicability 
limited to the cases where the accuracy required is not so 
much lower than just a few percent. Therefore, in the cal- 
culation of the geometrical and time acceptances, meth- 
ods which can use directly the data are very appreciated 
|13|14j . Obviously, such data based methods will face us 
with subtle effects, since the sample from which we are try- 
ing to extract the exposure, contains, in principle, the sky 
anisotropy which we wish to later estimate. The usual pro- 
cedure applied to build the exposure map from the data 
itself is to perform some kind of scrambling into the local 
angles and arrival times in order to wash out the possible 
sky anisotropics and retain only the local features [T5| . The 
role played by the time behaviour of the detector is there- 
fore central in calculating the sky exposure map since, on 
the one hand, detector long term variability induces fake 
anisotropic patterns on the sky, while on the other hand, 
large scale anisotropy can distort the time distribution of 
events registered by a stably running detector. As an ex- 
ample of the former effect, we know that a solar modula- 
tion superimposed on the top of a seasonal envelope, both 
genuine local weather effects, will give rise to an apparent 
sidereal modulation. And for the latter, it is widely known 
that any modulation in right ascension will appear as a 
sidereal variation in the time distribution of events. 



In this article we will deal only with the time variation 
of the exposure as this is usually the most delicate issue. 
Our method can be applied to any kind of detector with 
nearly time independent aperture, a typical property of 
ground arrays where the detector has, in principle, 100% 
duty cycle. Through the text, just for illustration purposes, 
we are going to work with a cosmic ray detector located at 
the same position of the Auger Observatory [16] , but the 
results should be valid for a large variety of detectors. 

The paper is organised as follows: we present the ba- 
sic underlying assumption of Julian day x solar factoriza- 
tion of the detector rate in section [21 and in section [3] we 
show how to recover the time detector distribution from the 
dataset itself if the possible sky anisotropics, through their 
sidereal amplitudes, are no larger than just a few percent, 
therefore avoiding the use of Monte Carlo simulation. The 
residual systematic effects in the coverage due to possible 
non-factorizable components of the detector rate is esti- 
mated in section |4] by using a model for the detector rate 
based on weather monitoring data taken at the Auger site. 
The robustness of coverage maps built under the factor- 
ization hypothesis is treated in section [5] by reconstructing 
some anisotropic patterns on the sky. We finally conclude 
in section [6l 

2. The Factorizable Acceptance Model (FAM) 

Let A be the detection efficiency of a surface array. This 
function depends on the time t, on the shower horizontal 
coordinates 6 and (p ^nd on a set of parameters, like the 
energy, which we will refer to as X. By writing the time 
as a function of two distinct variables: t = t{j, s), where j 
is the Julian day and s is the solar time, we can therefore 
represent the array efficiency as A{j, s, 9, 0, X). 

In the most general case, the efficiency will have an in- 
tricate form with non-trivial correlations among all the pa- 
rameters. For example, at a fixed energy, the weather in- 
duced time modulations might have a relative amplitude 
which depends on the zenith angle. In the following, one as- 
sumes that the domain spanned by the variables {0, 0, X} 
can be broken up into complementary regions over which 
A can be factorized as 

AU,s,9,<l>,X) = Atij,s) X An{9,(l>,X), (1) 

where At{j, s) and Anid, 4>, X) are both dimensionless, and 
the results are given for only one of such regions. In practice, 
to get the full picture, one does the same analysis over all 
regions required to be defined. 

In the factorizable acceptance model (FAM), we assume 
that the time variation of the efficiency is the product of a 
daily variation in solar time ?7d(s) and a seasonal variation 
in Julian days rjYij), so that one can write 

At{.j,s) ^7]y{j) xt]d{s). (2) 

The hypothesis behind this factorization is that what 
happens during a day is not exclusive to that particular 



day, but mainly correlated to the sun altitude in the sky 
(given by the solar time) in that day. The validity of such 
a hypothesis can be checked directly from the data set as 
will be shown in section [H 

3. Extracting the acceptance under the FAM 
hypothesis from the data set itself 

The goal of this section is to show how one can retrieve 
the time modulations of the detector acceptance from the 
data set itself, if one uses an integer number of data taking 
years and if a possible anisotropic pattern on the sky has 
a small amplitude, that is, a sidereal time modulation no 
larger than a few percent. For ultra high energy cosmic 
rays, this is in fact a pretty good assumptiorPl. 

Let J^{a, 6, X) be the cosmic ray flux from the direction 
(a, S) in equatorial coordinates, where a is the right ascen- 
sion and S is the declination, reaching the Earth per unit 
time, area, solid angle and per units of X. The variables a 
and S can be written in terms of the local horizontal coor- 
dinates 9, (f) and the time of detection t{j, s). Therefore, the 
following integral over one of the regions where the factor- 
ization ^ is valid 



A, 



6, X.)T{a{t, 0, (l)),S{t, 0, (/>), X)dXdn, (3) 



X n 



has dimension of number of particles per unit time and area. 
The rate of cosmic rays detected (over the region where 
is defined) as a function of time is then given by 



n{t) = n{j, s) = At{j, s) X $(0 X s.ffit), 



(4) 



where Seff{t) represents here the instantaneous effective 
area of the detector, since this might not be constant in 
time, as is the case of detectors which take data as their 
sizes are still growing by the installation of new detection 
units, or experimental problems which might black part of 
the detector, reducing its effective area. These effects can 
usually be corrected for and we will not be concerned with 
them in this paper. 

The function (f>(t) gives us the time dependence of the 
sky anisotropy and can be written in terms of the sidereal 
time tsid as 



(5) 



where e(tsid) is a function of period Dsid (a sidereal day) 
with null integral over this time interval and |e| << 1 and 
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^(tsid)dtstd, 



(6) 



^ We are assuming here as ultra high energy cosmic rays those with 
energies above ~ 10^^ eV. Therefore, the statement is based on 
the isotropy at PeV energies observed by KASCADE 2], the Auger 
results [3], and it remains valid even if wc take the AGASA dipolc 
[T] as real. 



Since a Julian day is about only 4 minutes longer than a 
sidereal day (a correction factor of about 1 -I- 1/365.25) 
and a calendar year being a little longer than the time it 
takes for t^id to span the whole sky in right ascension (same 
correction) , we have two useful approximations 

f^{tsid)dtsid ^ 0, and J e{ts^d)dtsid ^ 0, (7) 

where Ygoi and Dsoi are the durations of a solar year and a 
solar day, respectively. 

By explicitly writing the average values of ??y(j) and 
77£)(s), we have 

riY{j)=riY{l + fvij)) and ?7d(s) = 77^ (1 + //^(s)) , (8) 
and by definition the fractional amplitudes satisfy 

fvim^O. and j fDis)ds^O, (9) 

With all that in mind, we can write the rate of cosmic 
rays, corrected by the effective array surface as 



^ (j, s) 



: 77^.77^*0 [1 + e(i..d)] [1 + fy + Id + ,(10) 



where we have omitted the dependence of fy and fjj on j 
and s, respectively. Integrating over a solar day or a solar 
year and using eqs. ([7]) and (O, we obtain the functions 
T^'yU) = I 7^'(j,s)ds and 7^^(s) = / n'{j,s)dj 
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f\tsid)fD{s)ds 
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77- / (■{tsid)fY{j)dj 
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Y,al 



Vd{s).{12) 



We are specifically interested in the estimation of the cov- 
erage map for cosmic rays at energies above around 10^ eV, 
detected by cosmic rays detectors such as large surface ar- 
rays. As already mentioned, such arrays are quite stable in 
time, with typical fractional amplitudes for the time vari- 
ations in the detection efficiency which are no larger than 
10% (a total 20% variation) along a day or a year (see sec- 
tion|4]). Moreover, the sidereal time modulation at these en- 
ergies will not exceed the 5% level, a somewhat exaggerated 
scenario (see footnote on page [3]). Under these hypotheses 
the integral inside square brackets can be neglected, since 
it will be certainly below the 0.5% level. 

Therefore, one can write 

7^y(i) ~ ^QT]jjDsoiVYij), n'ois) ~ ^oVvYsoiVDis), (13) 

which contain only the detector efficiency functions, re- 
gardless of any anisotropy on the sky. The left hand side of 
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these two equations arc distributions which can be easily 
obtained from the data set itself as the histograms of phys- 
ical events as a function of the Julian day or the solar time. 
It is therefore possible to extract the detector information 
from the data, even in the presence of anisotropy. The ap- 
proximations in eq. (|13p are good up to order |e| |/y(_D) | (the 
product of the absolute values of e and fY{D))- For all those 
purposes where the required precision over the coverage es- 
timation lies below this level, the procedure outlined here 
does not apply. 

A coverage map can then be built from the data by us- 
ing the local reconstructed angles 9 and (j) and sorting the 
arrival times according to the quantity 



(14) 



where the proportionality factor depends on <I>o, rjy and 
rjjj, being therefore unknown, but completely unimportant, 
since it is a global constant factor not depending on time. 

4. The residual systematic due to non-factorizable 
components in the acceptance 

The trigger rate of a surface array depends on the weather 
conditions, since the changes in the pressure P and air den- 
sity p affect the shower development through the atmo- 
sphere. An increase in air density will reduce the Moliere 
radius, whereas the average shower age at ground increases 
with the atmospheric pressure. As shown in [17| . the Auger 
surface detector trigger rate can be fairly well modelled 
through a linear dependence with P and p. A fit to the data 
collected by the observatory spanning a 2-year period have 
produced: 

TZ{t) = i?o [1 + ap{P - Po) + Upipd - Po) + hp{p - pd)] , (15) 

where po = 1-055 kg m"'^ is the average reference density, pd 
is the average daily density, ap = (9±5) x 10~^ hPa"-'^, Op = 
(-2.68 ± 0.07) kg^hv? and bp = (-0.85 ± 0.07) kg-^m^. 

From now on, we will assume that the trigger rate of a 
surface array is well modelled by eq. (|15p. representing ap- 
propriately both its long term seasonal and daily modula- 
tions, as well as its short time fluctuations due to the un- 
derlying variations in pressure and densitjO]- Everywhere 
in the paper, we assume that the array is located at the 
same geographical position as the Auger detector (35.25° 

5, 69.25° W). Figure[I]shows the rate predicted by eq. p5)) . 
using the same weather monitoring data of reference [17] , as 
a function of the Julian day and the solar time. Two years 
are shown (2005-2006) and one can see the '^6% seasonal 
modulation (top) and daily (bottom) variation of ~2%. 

In order to estimate the possible residual systematic 
due to the presence of components in the trigger efficiency 




Fig. 1. Seasonal (top) and daily (bottom) rate variations according to 
the model represented by equation I I15I I shown in ref. 1171 to describe 
fairly well the rates observed in the Auger data set. The rate is 
assumed to depend linearly on the ground pressure and air density 
(see coefficients in the text and 1171 ). The latter is obtained from the 
temperature. Pressure and temperature measurements were taken 
from a weather station on site and two years of data arc shown: 
2005-2006. 

whose time behaviour cannot be factorized into a seasonal 
and a diurnal effect, we simulated an isotropic cosmic ray 
distribution, convolving it with a detector time efficiency 
given by eq. (jl5p and a simple geometrical efficiency given 
by sin 6 cos 6, which takes into account both the solid angle 
and the effective detector transversal area seen by a shower 
at zenitlO] angle 9. All the simulations in this paper include 
only showers with 6' < 60 degrees, consistent with the data 
sample used in [T7] . In the absence of anisotropic patterns 
from the sky, we know that all distortions observed in the 
arrival time distribution of the showers are due only to 
local effects. Therefore, the true coverage map can be built 
directly from the data set by performing a scrambling of 
the arrival times. The exposure sky map built in this way 
will be denominated here UTC, since the arrival times are 
usually given as UTC time. 

The validity of the factorization hypothesis can then be 
tested by comparing the UTC coverage with the one ob- 
tained by replacing the original arrival times by a value 
drawn according to eq. (fT4)) . The latter map will be called 
the FAM coverage. The coverage ratio FAM/UTC is shown 
in figure O binned in right ascension and declination. To 
beat the statistical fluctuations down and have access to 



Here, we assume the air density is completely determined by P 
and T according to p = p/RdryT, where Rdry = 287.05 J/kg/K is 
the specific gas constant for dry air, that is, we neglect the effect of 
humidity. 



For a real surface array, there are several additional effects which 
contribute to the coverage estimation. However, as already stated 
before, we are only interested here in the role of the time behaviour 
of the trigger efficiency. 
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Fig. 2. Ratio between the FAM and the UTC sky coverage map as 
a function of right ascension (top) and dechnation (bottom). The 
exposure maps were estimated by using simulated data with events 
following the pressure and temperature modulated rate modelled by 
eq. Hist . The error bars represent the statistical error only. Here, 
since there is no anisotropic pattern from the sky, the true coverage 
is represented by the UTC sky map. 

small systematic effects, these plots are based on large 
Monte Carlo samples (8 x 10^ events). One can see that 
the residual modulations are below the 0.2% level and arc 
still consistent with the statistical ffuctuations. Therefore, 
the coverage map of a surface array whose time trigger ef- 
ficiency can be fairly well described by a linear dependence 
on air pressure and density as given by cq. (I15p . can be ac- 
curately (that is, up to 0.2% level) built under the hypoth- 
esis of seasonal x diurnal effects factorization. 

5. The method performance under large scale 
anisotropy patterns 

Having validated the hypothesis of factorization in the 
last section, one can study the performance of the FAM 
and UTC coverages when in addition to local effects, the 
arrival time distribution is further distorted by large scale 
anisotropy patterns from the sky. 

A coverage map constructed from the data set itself 
which keeps the original arrival time distribution, like 
the UTC coverage, which essentially makes a scrambling 
of the arrival times, will always underestimate the am- 
plitude of a right ascension modulation, since all of the 
variation induced into the time distribution is absorbed 
into the coverage. The usual excess maps made by taking 
the signal/coverage ratio will therefore underestimate the 
modulation in right ascension. 

Figure [3] shows the time distribution used to build the 
coverage under the FAM hypothesis and the usual scram- 
bling of UTC shower arrival time, when an equatorial dipole 
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Fig. 3. Sidereal modulation for the UTC coverage (red squares) 
and the sidereal (black circles) and anti-sidereal (blue triangles) 
sidebands for the FAM coverage. The equal amplitudes sidereal and 
anti-sidereal sidebands for the FAM coverage are the results of the 
seasonal modulated solar variation. The dotted lines are fits to the 
respective sidebands with amplitudes: (0.10 ±0.07)% (sidereal) and 
(0.11 ±0.07)% (anti-sidereal). We have artificially introduced a shift 
in the phase of the anti-sidereal band for better visualisation, since 
this one is practically in phase with the sidereal band. 

of amplitude 2% is present. The dipolar pattern induces a 
large sidereal contamination into the UTC coverage as can 
be seen in figure [3l In fact, this is one of the most unfavor- 
able cases for a UTC like coverage map, since it is a pure 
right ascension distortion, creating a large modulation in 
the events arrival time distribution. The concept of anti- 
sidereal time was introduced some time ago by Farley and 
Storey [H] in order to study long term behaviour of cosmic 
ray detectors. The authors made use of the fact that a sea- 
sonal modulated diurnal variation (both assumed to be har- 
monic) could be obtained by the interference of sidebands 
of equal amplitudes and slightly different periods, in which 
the slightly higher frequency (with respect to a solar fre- 
quency) sideband is the well known sidereal band, whereas 
the slightly lower frequency sideband was called the anti- 
sidereal component^- Therefore, even pure local seasonal 
and diurnal modulations will give rise to an apparent side- 
real anisotropy, but these can then be identified by looking 
for the associated equal amplitude anti-sidereal sideband. 
The time distribution of events used in the FAM coverage 
shows in fact such twin sidebands with fitted amplitudes 
of about 0.1%. The amplitude of the sidebands is given ap- 
proximately by the product of the amplitudes of the sea- 
sonal and the diurnal components. The fitted amplitudes 
to the rates of figurc[l]imply an expected 6% x 2% = 0.12% 
sidereal (anti-sidcrcal) modulation, which is in pretty good 



* An anti-sidereal year has approximately one day more than a solar 
year and a sidereal year has one day less. 
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Fig. 4. Distribution of reconstructed amplitudes of an equatorial 2% 
amplitude dipole via the method of reconstruction of pseudo-multi- 
polar coefficients for non- uniform and incomplete sky coverage 1221 . 
Two different coverages: UTC (shaded blue) and FAM, and used as 
input to the reconstruction are compared. The central values and 
the standard deviations for the fitted Gaussians are shown on the 
top of each histogram. Each one of the fOOO samples shown have 
5 X 10^ events. 

agreement with figure [3] Therefore, by assuming a Julian 
day X solar time factorization of the rate, essentially none 
of the right ascension modulation leaks into the final ex- 
posure sky map. It is worth stressing, though, that such 
a cancellation of sidereal contribution from sky anisotropy 
can only be achieved by using an integer number of years 
as already discussed in section[3l We stress as well that un- 
like Farley and Storey which assume a harmonic form for 
the seasonal and diurnal modulations, in the FAM method, 
their analytical forms are extracted directly from the data. 

A more quantitative idea of the bias introduced into the 
coverage by a simple scrambling of the original arrival times 
is given in figured The plot shows the distribution of recon- 
structed amplitudes for 1000 realisations of the 2% equa- 
torial dipole. The reconstruction of dipolar patterns with 
a full sky observatory was first analysed in |19| with later 
generalisations to the case of partial sky coverage presented 
in [20] and [21]. In figurelH the reconstruction is performed 
by retrieving the set of multipolar coefficients of the events 
directions map for a non-uniform and, in this case, even in- 
complete sky exposure. The central idea of the method is 
that the pseudo multipolar coefficients obtained with non- 
uniform partial sky coverage are related to the true ones by 
a convolution operation whose kernel is a function of the 
detector window function (the coverage map) P^j . 

From figurelUwe see a clear underestimation of the dipole 
amplitude when the reconstruction is done by using the 
UTC sky map, with the average value (1.28 ± 0.35)% being 
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Fig. 5. Dipole amplitudes reconstructed through the dcconvolution 
of the multipolar coefficients for non-uniform and incomplete sky 
coverage 1221 as a function of the declination of the dipole axis. 
FAM (black circles) and UTC (red squares) based coverage maps are 
compared. The true 2% amplitude is represented by the horizontal 
dashed line. 

2(7 away from the true value. With the FAM coverage, in 
turn, one is able to retrieve the amplitude within la. The 
relative statistical uncertainty of 11% in the amplitude for 
the FAM case is also smaller than the 27% corresponding 
to the UTC coverage. 

It is interesting to compare these coverage based recon- 
structions with a coverage independent one. For a pure right 
ascension (that is sidereal) modulation, it is straightfor- 
ward to estimate the sidereal amplitude through the power 
distribution at the sidereal frequency of the arrival time 
Fourier transform as described in [12j . The reconstruction 
is done in a two step way. Firstly, the normalised Fourier 
coefficients of a set of N modified times are calcu- 

lated at frequency / (tj = 27r/) as a(/) = c(/) + is{f) 



s{f) 



W 



Wi siniujti) (16) 



with W = (^i = l)i where the sums run all over 

the N events. An estimation for the annual (diurnal) am- 
plitude pj^(d) = ^s2(/y(d)) + c^{fy(d)) and phase (j)y(d) = 
atan(s(/j,((;)), c{fy(^cL))) is then obtained. Finally, the Fourier 
amplitudes in (I16p are recalculated, now with weights 

= [(1 +pyCos{ujU - (l>y)) X (1 +pdCos{ujU - M)]^^ (17) 

to get the final sidereal amplitude and phase. Such a Fourier 
analysis, totally independent of the coverage map, gives an 
amplitude of (1.6±0.2)% for the same 2% equatorial dipole, 
a result somewhere in between the reconstructions through 
the multipolar coefficients. 

A more systematic test of the performance of the FAM 
and UTC based coverage maps is shown in figure [S] where 
the reconstructed amplitudes and their associated errors 
are shown for dipoles with axis at different declinations 
(and fixed amplitude of 2%). By varying the angle between 



^ It was shown in I12| that by correcting the arrival times by the 
right ascension phase of the event with respect to the local sidereal 
time, one is able to observe modulations on smaller scales and resolve, 
for example, the sidereal and diurnal frequencies as long as data is 
gathered for a sufficiently long period of time. All the simulations in 
this section are for a 5-year period. 



its axis and the equatorial plane, one changes the ampli- 
tude of the induced sidereal variations. One can see that, in 
average, the amplitude is more accurately retrieved when 
the FAM coverage map is used as input to the multipolar 
reconstruction. In the whole range between -90° and 90°, 
the UTC coverage leads to a mean amplitude 2.9cr away 
from the true value, whereas a 1.3a estimation is achieved 
with the FAM sky map. 

6. Conclusions 

We have introduced a new method for estimating the ex- 
posure sky map of a cosmic ray detector, based on the as- 
sumption of factorization of the time behaviour of such a 
detector into seasonal (represented by the Julian day) and 
diurnal (represented by the solar time) variations. Such a 
hypothesis is motivated by the current knowledge on the 
modulations of the detection rate of large surface arrays 
caused by weather effects, such as pressure and air temper- 
ature, which affects both the longitudinal and lateral devel- 
opments of particle cascades through the atmosphere. By 
using an integer number of data taking years, one is able 
to cancel the sidereal contribution of sky anisotropy to the 
Julian day and solar time histrograms, allowing to extract 
the detector time behaviour from the data set itself. 

In |17j . a model was shown to describe fairly well the 
yearly and diurnal rate modulations observed in the Auger 
data set, where the rate is assumed to vary linearly with 
pressure and air density, with the proportionality coeffi- 
cients being fitted directly to the number of showers de- 
tected by the surface array and to the atmospheric mon- 
itoring data taken regularly on site. Using this model, we 
have shown that possible residual systematic effects caused 
by the presence of non-factorizablc terms into the long term 
detector behaviour will induce fake large scale anisotropy 
modulations on the sky no larger than 0.2%, being there- 
fore negligible for a typical 1% desired accuracy into the 
coverage estimation. 

Coverage maps built under the factorization hypothesis 
were shown to be more suitable when reconstructing dipo- 
lar anisotropic patterns through the recovery of its multi- 
polar coefficients, providing amplitudes with a mean sys- 
tematic shift of 1.3(7 with respect to the true value as one 
changes the declination of the dipole axis. The amplitudes 
reconstructed from a coverage where one simply scrambles 
the showers arrival times are strongly biased, with values 
tipically around 3a away from the true values. A cover- 
age independent estimation of the equatorial dipole, such 
as the one provided by the Fourier analysis of the arrival 
times and detector running for a 5-year period gives an am- 
plitude somewhere in between the multipolar method with 
UTC and FAM input coverages, with a systematic under- 
estimation of the amplitude at the 2a level. 

There are some similarities behind the factorization hy- 
pothesis with the analysis performed in |18j, but whereas 
Farley and Storey assume a particular harmonic depen- 



dence for both the seasonal and diurnal modulations in or- 
der to fit their respective amplitudes, the time distribution 
used to build the FAM coverage is taken from the data set 
itself. 
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